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Abstract 

Although two-dimensional methods provide accurate 
predictions of contact stresses and bolt load distribution 
in bolted composite joints with multiple bolts, they fail 
to capture the effect of thickness on the strength 
prediction. Typically, the plies close to the interface of 
laminates are expected to be the most highly loaded, 
due to bolt deformation, and they are usually the first to 
fail. This study presents an analysis method to account 
for the variation of stresses in the thickness direction by 
augmenting a two-dimensional analysis with a one- 
dimensional through the thickness analysis. The two- 
dimensional in-plane solution method based on the 
combined complex potential and variational formula- 
tion satisfies the equilibrium equations exactly, and 
satisfies the boundary conditions and constraints by 
minimizing the total potential. Under general loading 
conditions, this method addresses multiple bolt configu- 
rations without requiring symmetry conditions while 
accounting for the contact phenomenon and the inter- 
action among the bolts explicitly. The through-the- 
thickness analysis is based on the model utilizing a 
beam on an elastic foundation. The bolt, represented as 
a short beam while accounting for bending and shear 
deformations, rests on springs, where the spring 
coefficients represent the resistance of the composite 
laminate to bolt deformation. The combined in-plane 
and through-the-thickness analysis produces the bolt/ 
hole displacement in the thickness direction, as well as 
the stress state in each ply. The initial ply failure 
predicted by applying the average stress criterion is 
followed by a simple progressive failure. Application of 
the model is demonstrated by considering single- and 
double-lap joints of metal plates bolted to composite 
laminates. 
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Introduction 

Bolts provide the primary means of connecting com- 
posite parts in the construction of aircraft and aerospace 
vehicles. The main disadvantage of bolted joints is the 
formation of high stress concentration zones at the loca- 
tions of bolt holes, which might lead to a premature 
failure of the joint due to net-section, shear-out, or bear- 
ing failures, or their combinations. The stress state in a 
bolted joint is dependent on the loading conditions, 
dimensions, laminate stacking sequence, bolt clamp-up 
forces, bolt location, bolt flexibility, bolt size, and bolt- 
hole clearance (or interference). A substantial number 
of experimental, analytical, and numerical investiga- 
tions have been conducted on the stress analysis of 
bolted laminates. The study by Kradinov et al. 1 pro- 
vides an extensive and detailed discussion of earlier 
investigations. In order to eliminate the shortcomings of 
the previous analyses, Kradinov et al. introduced a two- 
dimensional numerical/analytical method to determine 
the bolt load distribution in bolted single- and double- 
lap composite joints utilizing the complex potential- 
variational formulation. This method addresses multiple 
bolt configurations without requiring symmetry condi- 
tions while accounting for the contact phenomenon and 
the interaction among the bolts explicitly under bearing 
and by-pass loading. The contact stresses and contact 
regions are determined through an iterative procedure 
as part of the solution method. 

Although this two-dimensional approach provides an 
accurate prediction of the contact stresses and bolt load 
distribution, it fails to capture the effect of thickness on 
the failure prediction. In addition to the head and nut 
shapes and the applied bolt torque, the stacking se- 
quence considerably influences the stress state in each 
ply of the laminate. Thus, an adequate representation of 
the ply load variation through the thickness is critical 
for the failure prediction of composite laminates at the 
bolt-hole boundary. 

This study presents an analysis method to account 
for the variation of stresses in the thickness direction by 
augmenting the two-dimensional analysis by Kradinov 
et al. 1 with a model of a bolt on an elastic foundation, as 
suggested by Ramkumar et al. 2 The bolt, represented as 
a beam while accounting for bending and shear defor- 
mations, rests on springs, where the spring coefficients 
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represent the resistance of the composite laminate to 
bolt deformations. The values of the spring coefficients 
depend on the fiber orientation of the laminate plies; for 
isotropic plates, the spring coefficients are defined by a 
constant value. The present analysis produces the bolt/ 
hole displacement in the thickness direction and the 
stress state in each ply. Failure load and associated fail- 
ure modes of net-section, bearing, and shear-out for 
composite bolted joints are predicted based on the aver- 
age stress criterion of Whitney and Nuismer 3 for first 
ply failure, followed by a simple progressive failure cri- 
terion as suggested by Ramkumar et al. 2 

The applicability of this method is demonstrated by 
considering single- and double-lap joints of laminates 
with a varying number of bolts. In addition to the deter- 
mination of the contact stresses and the bolt load distri- 
butions, the failure load is investigated by applying a 
progressive failure procedure based on the average 
stress failure criterion. 

Problem Statement 

The geometry of bolted single- and double-lap joints 
of composite laminates is described in Fig. 1. Each joint 
can be subjected to a combination of bearing, by-pass, 
and shear loads. Each laminate of the single- and 
double-lap joints, joined with L number of bolts, can be 
subjected to tractions and displacement constraints 
along its external boundary. The thickness of the 
laminates is denoted by h k . As illustrated in Fig. 2, the 
hole radius in the k th laminate associated with the I th 
bolt, a k t (which is slightly larger than the bolt radius, 
R f ), leads to a clearance of 8 k i . The ranges of the 
subscripts are specified by and £-\ ,...,£, 

with K and L being the total number of laminates and 
bolts, respectively. The bolt radius remains the same in 
each laminate; however, the radii of the holes associ- 
ated with the same bolt are not necessarily the same. 

The extent of the contact region is dependent on the 
bolt displacement deformation of the hole boundary, 
and the clearance. The presence of friction between the 
bolts and the laminates is disregarded. Each laminate 
with a symmetric lay-up of N k plies can have distinct 
anisotropic material properties. Each bolt can also have 
a distinct stiffness, and the explicit expressions for bolt 
stiffness for a single- and double-lap joint, as well as 
the general lap configurations, are derived in Kradinov 
et al. 1 

The problem posed concerns the determination of the 
extent of the contact zones, the contact stresses and the 
bolt load distribution under general loading conditions, 
the bolt/hole deformation, and the stress state in each 
ply, and thus the joint strength. 





Fig. 1 Geometric description of single- and double- 
lap bolted joints. 



Fig. 2 Position of a bolt before and after the load is 
exerted. 


Solution Method 

In-Plane Analysis for Contact Stresses and 
Bolt Load Distribution 

The coupled complex potential and variational for- 
mulation introduced by Kradinov et al. 1,4 is employed 
to determine the two-dimensional stress and strain 
fields required for the computation of the contact 
stresses and contact regions, as well as the bolt load dis- 
tribution. This in-plane analysis is capable of accounti- 
ng for finite laminate planform dimensions, uniform 
and variable laminate thickness, laminate lay-up, inter- 
action among bolts, bolt torque, bolt flexibility, bolt 
size, bolt-hole clearance and interference, insert dimen- 
sions, and insert material properties. Unlike the finite 


2 

American Institute of Aeronautics and Astronautics 


element method, it alleviates the extensive and expen- 
sive computations arising from the non-linear nature of 
the contact phenomenon. Also, the method is more suit- 
able for parametric study and design optimization. 

Although this two-dimensional analysis provides 
accurate in-plane stresses in each laminate and bolt load 
distribution, it assumes no variation of stresses through 
the laminate thickness. This assumption might lead to 
erroneous results in the strength prediction of bolted 
joints because of the pronounced influence of through- 
the-thickness stress variation at the bolt location as 
discussed by Ramkumar et al. 2 

Through-Thickness Analysis for 
Bolt/Hole Deformation 

In conjunction with a two-dimensional in-plane bolt- 
ed joint analysis, Ramkumar et al. 2 suggested a model 
utilizing a beam on an elastic foundation in order to 
include the variation of stresses in the thickness direc- 
tion of the bolted joint. The bolt rests on springs, where 
the spring constants represent the resistance of the lami- 
nate to bolt deformation. The spring constants corre- 
spond to the modulus of each ply through the thickness 
of the hole boundary. Their values depend on the ply 
orientation of the laminate. For isotropic plates, the 
spring constants have a uniform value. 

As the bolt bends, the plies are loaded differently 
near the hole boundary based on their orientation and 
location. As shown in Fig. 3, the plies close to the inter- 
face of adjacent laminates exhibit significant deforma- 
tion. As shown in Fig. 4, the beam representing the bolt 
rests on an elastic foundation whose modulus is 
represented by the stiffness of each ply, k k J , in the 
laminate. The superscript f and subscripts k and i 
denote the specific bolt, the laminate, and the ply num- 
bers, respectively. Also, the bolt is subjected to con- 
straints at the head and nut locations through rotational 
stiffness constants, k[ l) and k [ l) , in order to include the 
effect of head and nut shapes and bolt torque. 

Free-body diagrams of the laminates at the I th bolt 
in a single- and double-lap joint and the end conditions 



Fig. 3 Typical bolt deformation in 


and slope continuity conditions in the presence of both 
bending and shear deformations are shown in Fig. 5. In 
accordance with the typical bolt deformation illustrated 
in Fig. 3, the force exerted by the I th bolt on the k th 
laminate, Pjf ] (obtained from the two-dimensional in- 
plane analysis), is enforced as a shear force, V k (£) , at the 
interface of the adjacent laminates. At the interface, the 
continuity of the bending slopes, x// k } = y/ k \ , is also 
enforced while permitting the laminates to displace. At 
the head and nut locations of the bolt, the shear force 
values are set to zero, and the rotations (slopes) are 
dictated by rotational stiffness constants, k { h l) and k [ l) , 
depending on the bolt type, presence of washers, and 
the applied bolt torque. 

Finite Element Analysis . The bolt/hole displace- 
ments through the laminate thickness are obtained by 
discretizing the bolt with beam elements that account 
for bending and shear deformations. The bolt discreti- 
zation is based on the discrete nature of the ply stacking 
sequence. Along its thickness, each ply is discretized 
with two beam elements. For both single- and double- 
lap joints, the number of elements and the number of 
nodes in relation to the number of plies in the laminate 
are described in Fig. 6. In the discretization process, a 
node located in the middle of each ply is attached to a 
spring element representing the ply stiffness (Fig. 6). 

The derivation of the stiffness matrix composed of a 
two-noded beam element (Timoshenko’s zeroth-order 
shear deformable beam theory) and a linear spring ele- 
ment is presented in the Appendix. 

Associated with the k th laminate and the I th bolt, 
each node is assigned a deflection, A k \ = A[ n (z . ) , and 
a rotation, (j) k \ - (j) k ] (z . ) , with the subscript j repre- 
senting the node number. In the finite element 
formulation, the rotations of the internal nodes are 
statically condensed in terms of the nodal displace- 
ments and the end node rotations. The positive direc- 
tions of the deflections and rotations are shown in Fig. 
7. The details of the condensation procedure are also 
explained in the Appendix. 



(a) single- and (b) double-lap joint. 
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Fig. 4 Bolt on an elastic foundation model in a (a) 
single- and (b) double-lap joint. 
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Fig. 5 Free-body diagrams of laminates for a (a) 
single- and (b) double-lap joint. 
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Fig. 6 The finite element model of a bolt in a (a) single-and (b) double-lap joint. 
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Fig. 7 Bolt discretization in the k th plate after static condensation. 
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Spring Stiffness Coefficients . As suggested by 
Ramkumar et al. 2 , the translational spring stiffness 
coefficients, k k j , representing the i th ply of the k th 
laminate near the I th bolt are approximated by 


h ’ 1 rf 


i = \,N k (1) 


employed in the prediction of the initial ply failure load, 
F^j IN , and the corresponding failure mode. 

As shown in Fig. 2, the maximum hole enlargement, 
y ( k l) , is defined as the absolute value of the difference 
between the radial displacements, u[ l) (r = a k - 6 X ) 
and u k \r = a kl ,9 -0 2 ) ,of points 1 and 2 on the hole 
boundary in the direction of the bolt load 


where p k J is the load exerted by the I th bolt on the 
I th ply of the k th laminate and A 0 represents the maxi- 
mum hole enlargement of the I th hole in the k th lami- 
nate. The number of plies in the k lh laminate is denoted 

by N k . 

As part of the two-dimensional in-plane bolted joint 
analysis, the load exerted by the bolt on the i th ply of 
the k th laminate near the I th hole, p k J , is computed as 

P'S = psJM’J < 2 > 

where p[l] i)x and p{l] i)y represent its components in the 
x- andy-directions. These components are computed by 
integrating the radial stresses in each ply as 

pfk,i )x = a k ,t J AY ( r = ,0) cose de (3a) 

0 

and 

p{Z)y = a k,t J ( r = a k,i ,e) sine do (3b) 

0 

in which a ki is the radius of the I th hole in the k th 
laminate, and <y (k / ] (r = a kJ , 6) represents the radial 
stress distribution in the i th ply of the k th laminate near 
the I th hole. 

Under plane-stress assumptions, the stress and strain 
components are related by 

(4) 

in which Q k i represents the reduced stiffness matrix 
for the i th ply of the k th laminate. The Cartesian stress 
components in the i th ply of the k th laminate and the 
Cartesian strain components, uniform through the 
thickness of the k th laminate, are included in the vec- 
tors of a {k,l) and z {k) . This stress state in each ply is 


YT = | Ut\r = a kJ ,e = e 2 )-uf(r = a kt ,e = 6>,)| (5) 

where the radial displacements are obtained from the 
in-plane bolted joint analysis. The maximum hole- 
enlargement, y k ] , can be different for each laminate of 
the bolted joint, but is uniform through the thickness. 
Also, it is specific to each bolt-hole in the laminate 
because it is dependent on the deformation response 
and the bolt load distribution. The head and nut 
rotational stiffness coefficients, k[ l) and k[ l) , respec- 
tively, have values close to zero for free-end conditions 
and to infinity for protruding head bolts under high 
torque. The stiffness matrix becomes singular if these 
coefficients approach zero. 

The analysis results include displacements and rota- 
tions, A { k ] j and (j) k \ , at the j th node, as well as the 
spring forces, f k . , at the I th ply of the k th laminate 
near the I th hole. The effect of through-the-thickness 
variation is invoked in the in-plane stress analysis by 
considering the spring forces, f k . , as the corrected ply 
loads. 

Progressive Failure Prediction 

There are three major failure modes in bolted com- 
posite lap joints: net-section, shear-out, and bearing 
(Fig. 8). The net-section failure is associated with fiber 
and matrix tension failure and shear-out and bearing 
failures are associated with fiber and matrix shear and 
compression failures, respectively. Failure in bolted 
laminates can be predicted by evaluating either the spe- 
cific stress components or their interaction at charac- 
teristic distances from the hole boundary. Although any 
one of these criteria is applicable to the prediction of 
the failure of a laminate or a ply, values of the charac- 
teristic distances and the unnotched strength parameters 
of the material are scarce. 



Net-tension Failure Shear-out Failure Bearing Failure 

Fig. 8 Primary failure modes in bolted composite joints. 
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The point and average stress criteria introduced by 
Whitney and Nuismer 3 disregarded the interaction 
among the stress components. However, they are wide- 
ly used in engineering practice for predicting the failure 
stress and failure modes because of their well-estab- 
lished values of the characteristic distances. 2,5 " 6 Both of 
these criteria predict net-section, shear-out, and bearing 
failures when the stress components at specific loca- 
tions reach their corresponding unnotched strength 
levels. The characteristic distances of a™ for net- 
section, a b o r for bearing, and a° for shear-out failures, 
as well as the shear-out and net-section planes (denoted 
by the n and s lines), are shown in Fig. 9. According to 
the point stress criterion, the net-section failure occurs 
when the normal stress, cr ss , at a distance a n o s from the 
bolt-hole boundary along the net-section plane reaches 
the unnotched tensile strength of a ply, X t . If cr ss at a 
distance a b o r from the bolt hole boundary reaches the 
unnotched compressive strength of a ply, X c , bearing 
failure occurs. Shear-out failure occurs when the shear 
stress, <j ns , at a distance a o ° from the bolt-hole 
boundary along the shear-out planes reaches the un- 
notched shear strength of a ply, X s . The average stress 
failure criterion is based on the average values of the 
corresponding stress components over the characteristic 
distances of (0,a” s ), (0,^ r ),and (0, a s 0 ° ) . 

Under the specified external loading, the ratios 
Cj k J £) {j - ns, br, so) of the unnotched strength param- 
eters to the average stresses associated with the net- 
section (j = ns) , bearing (j = br) , and shear-out 
(j = so) failure modes at the i th ply of the k th 
laminate near the I th hole are defined in the form 

C (b, P = — t 7r — for net-section failure mode (6a) 

ns, i —(£,*) v ' 

^ ss 



Fig. 9 Characteristic distances for point and aver- 
age stress failure criteria. 


Clf - _ c (k /) for bearing failure mode (6b) 

^ ss 

Cpf = _ s for shear-out failure mode (6c) 

® ns 

in which cff s ,l) and <j {b,l) are averaged normal and 
shear stresses over the characteristic distances (0, a™ ) , 
(0, a br ) , and (0 ,<). 

The initial ply failure load, F k f IN , and its associated 
failure mode are established by 

FkT = min [ c u l) ) Pk,i ’ (j = ns ’ br > ' s ' 0 ) (7) 

After the initial failure, a ply is assumed to continue 
sustaining the applied load according to a bilinear 
behavior, shown in Fig. 10. The value of the ultimate 
ply failure load, F k p UL , is defined by 

f (() ul = HF (e)m ( 8 ) 

where the factor H varies as 1.02, 1.50, and 1.12 for 
net-section, bearing, and shear-out failure modes, 
respectively, as suggested by Ramkumar et al. 2 

Due to the bilinear ply load behavior, the applied 
joint load is increased incrementally while predicting 
ply failure subsequent to the initial ply failure. At each 
load increment, the corrected ply loads, f k p , are 
compared to the initial and ultimate ply failure loads of 
F k p w and F k p uL , which are predicted according to the 
average stress criterion of two-dimensional analysis. 
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For an undamaged ply, if the corrected ply load of 
fff exceeds the corresponding initial ply failure load 
of F^ m , the ply experiences initial failure. Accord- 
ingly, as suggested by Ramkumar et ah, 2 the initial ply 
stiffness of k^J is reduced to kfj . The reduced ply 
stiffness, kfj , is defined by 

%?i=<*k% (9) 

in which the parameter a is assumed to be 0.1. For a 
damaged ply, if the corrected ply load of fff exceeds 
the corresponding ultimate failure load of F^ UL , the 
ply experiences total failure. Consequently, the ply 
stiffness is reduced to zero. 

Based on the bilinear behavior of the ply load shown 
in Fig. 10, the ply load at the i th ply of the k th laminate 
near the I th bolt can be expressed in terms of the ply 
displacement as 

4?=*^ (10a) 

for an undamaged ply, as 

fS=(^-^y^ + k^ (iob) 

for a damaged ply, and as 

//,?= o (10c) 

for a totally damaged ply, where j denotes the node 

associated with the translational spring element repre- 
senting the i th ply. 

When a ply fails, the adjacent plies share the load 
released by the failed ply. Thus, the failure propagates 
from ply to ply until the total failure of the laminate. 
The ultimate joint failure load is defined as the joint 
load that results in the ultimate failure of half of the 
plies at a particular bolt location. The minimum of the 
failure loads predicted for each bolt establishes the 
strength of the joint. This type of progressive failure 
analysis can be employed in conjunction with any one 
of the available failure criteria. 

Numerical Results 

The capability of this combined in-plane and 
through-the-thickness bolted joint analysis is 
demonstrated by considering single- and double-lap 
bolted joints joining metal to composite laminates with 
one, three, and four bolts as shown in Figs. 1 1-13. 

The material properties, stacking sequence and thick- 
ness of the plates are the same as those considered by 



Fig. 11 One-bolt single-lap joint geometry and 
loading. 



Composite laminates 


Fig. 12 Three-bolt double-lap joint geometry and 
loading. 
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Ramkumar et al. 2 The metal plates are made of alumi- 
num with Young’s modulus is fl =10.1Msi and 
Poisson’s ratio v a = 0.3 . The bolts are of steel with a 
Young’s modulus is 5 =30.0Msi and Poisson’s ratio 
v s = 0.3 . Although not a limitation of the analysis 
method, in these configurations, the bolt and hole 
diameters are equal, leading to zero clearance. The 
aluminum plates have a thickness of 0.31 in. The thick- 
ness of the laminate is 0.12 in with stacking sequence 
of [(45°/0 o /-45 o /0 o ) 2 0 o /90°] 5 . The material prop- 
erties for each ply are specified as E L =18.5Msi , 
E t = 1.9Msi , G lt = 0.85Msi , and v LT =0.3 . The 
high value of torque applied on the protruding bolt-head 
is specified by the head and nut rotational stiffness 
coefficients of k h =10 12 lbs-in and k n =10 12 lbs-in. 
The failure prediction is performed by employing the 
average stress criterion along with a bilinear stiffness 
reduction after the initial failure of each ply. The char- 
acteristic length parameters for the average stress fail- 
ure criterion are taken as a” s = 0.1 in, < 2 ™ =0.08 in, and 
=0.025 in. The unnotched strength parameters of the 
ply for each orientation in the stacking sequence are 
given in Table 1. 


Table 1 Unnotched strength values in Y-direction. 


Ply 

orientation 

(degree) 

X t , Net- 
section 
tensile (ksi) 

X C’ 

Bearing 

(ksi) 

X s> 

Shear- 

out 

(ksi) 

0 

230.0 

320.0 

17.3 

45 

40.0 

56.0 

95.0 

-45 

40.0 

56.0 

95.0 

90 

9.5 

38.9 

17.3 


As part of the finite element modeling, the section of 
the bolt in contact with the composite laminate is dis- 
cretized with 41 nodes in order to represent 20 plies of 
the laminate lay-up. Because the aluminum plate is 
thicker than the laminate, it is discretized with 8 1 nodes 
leading to 40 layers of aluminum. 

One-Bolt Single-Lap Metal to Composite Joint 

The geometrical parameters shown in Fig. 11 are 
defined by W =1.875 in, d =0.3125 in, L x =3.6 in, 
L 2 =4.4375 in, and s =0.9375 in. The initial applied 
load of P = 1875 lbs is uniformly distributed along one 
edge of the aluminum plate while the other end of the 
laminate is constrained. The variations of radial and 
tangential stresses around the hole boundary in alumi- 
num and composite plates are shown in Figs. 14 and 15, 
respectively. These figures demonstrate the capability 
of the two-dimensional analysis to capture the stress 
concentrations and provide the contact region around 



Fig. 14 Stress variation around the hole boundary 
in an aluminum plate. 



Fig. 15 Stress variation around the hole boundary 
in a composite laminate. 

the bolt hole. The segment of the radial stresses with 
negative values establishes the contact region between 
the bolt and hole boundary. As expected, there are zero 
shear stresses on the hole boundary because of the 
absence of friction. 

Based on the in-plane stress analysis, the maximum 
hole enlargements are computed as y x ] = 6.808x1 O' 4 in 
and y 2 l) =2.891xl0' 4 in for the aluminum and lami- 
nate, respectively. Invoking these values in Eq. (1), the 
stiffness of the spring representing the aluminum layer 
has a value of k[ l) t =69,597 lb/in with i = 1, 40 . The 
spring stiffness value for each ply of the laminate is 
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calculated as kf m = 34, 286 lb/in, = 52,586 lb/in, 

=34,654 lb/in’ and k® = 8, 446 lb/in, where the 
subscripts m, and g represent 45°, 0°, -45°, and 90° 
plies, respectively. 

The variation of the nodal displacements, Afj with 
(y' = l,81) and A with (7 =1,41) , illustrates the 
bolt/hole deformations in Fig. 16. As observed in this 
figure, the deformations in the composite laminate are 
larger than those in the metal plate as dictated by the 
material properties and laminate thickness. As expect- 
ed, the specified large values for head and nut rotational 
stiffness coefficients, k h and k n , result in zero slopes at 
the ends of the bolt. The maximum bolt/hole deforma- 
tions occur at the interface of the two plates, indicating 
the location of the major load transfer, as reflected in 
Fig. 17, which depicts the variation of the load distri- 
bution through the thickness of the joint. As expected, 
the load distribution through the thickness of aluminum 
plate varies continuously. However, the ply loads corre- 
sponding to the composite laminate change abruptly, 
depending on the fiber orientation. This behavior is dic- 
tated by the material property discontinuity in the thic- 
kness direction resulting in a different stress state in 
each ply. 

As presented in Table 2, the initial ply failure is 
predicted at a load level of P = 3,656 lbs, with a net- 
section failure mode in ply number 10 with a 90° fiber 
orientation. As the applied joint load is increased incre- 
mentally, the plies with a 90° fiber orientation continue 
failing in the net-section failure mode. Their failure is 
followed by a mixture of ±45° and 0° plies in the net- 
section and bearing failure modes, respectively. The 
load increments resulting in no failure have been 
omitted in Table 2. 

At load increment 56, ply number 1 with a 45° fiber 
orientation ultimately fails at a load level of P- 
5,336 lbs. This ply failure is followed by eleven differ- 
ent ply failures at the same load level. Therefore, the 
ultimate joint failure is reached at load increment 67 at 
a load level of 5,336 lbs. This prediction is in accep- 
table agreement with the experimental measurement of 
4,910 lbs reported by Ramkumar et al. 2 

Three-Bolt Double-Lap Metal to Composite Joint 

The geometrical parameters shown in Fig. 12 are 
defined by W - 2.05 in, L x = 3.6 in, L 2 = 4.525 in, v 1 
= 1.025 in, s 2 —1.0 in, s 3 = 0.9 in, s 4 =1.8in, \- 
1.025 in, h 2 - 0.5 in, h 2 - 0.4 in, and d - 0.3125 in. 
The initial joint load of P = 205 lbs is applied to the 
aluminum plate while the ends of the composite 
laminates are constrained. The maximum hole enlarge- 
ment values associated with each bolt hole are com- 
puted from the two-dimensional analysis and are pre- 


sented in Table 3, and the spring stiffness values for 
each ply are in Table 4. 

The through-the-thickness variation of the ply loads 
near bolt number 3 is shown in Fig. 18. The corre- 
sponding bolt/hole deformations are depicted in Fig. 19. 
As observed in these figures, the most pronounced 
deformation occurs in plies located along the plate 
interfaces. Both deformations and ply load distributions 
are identical for composite laminates due to the pres- 
ence of symmetry in the material and geometry. 

Bolts 2 and 3 exert higher loads on the composite 
than bolt 1 . The sequence ply failure loads and modes 
associated with each bolt are different because of the 
different strain states in the laminate near each bolt 
hole. 

As presented in Table 5, the initial ply failure near 
bolt 1 occurs at a load level of 22,402 lbs, in ply num- 
ber 19 with a 45° fiber orientation, in the shear-out 
failure mode. Part of the laminate near bolt 1 becomes 
unstable at load increment 23, corresponding to a load 
of 24,745 lbs, in ply 10 with a 90° fiber orientation, in 
the net-section ultimate failure. At this load level, 
seventeen more failures occur in the composite lami- 
nate before the laminate is assumed to ultimately fail at 
load increment 39. 

As presented in Tables 6 and 7, the initial ply failures 
near bolts 2 and 3 occur at 13,124 lbs and 12,769 lbs, 
respectively, in ply 1 1 with a 90° fiber orientation, in 
the net-section failure mode. The progress of failure 
near bolt 2 is presented in Table 6. Starting at load 
increment 46 and until 64, failure occurs for nineteen 
increments in different plies at a load of 18,225 lbs, and 
the joint can still carry more load. Finally, ultimate 
failure of the joint occurs at load increment 69, corre- 
sponding to a load level of 19,155 lbs, in ply 2 with a 0° 
fiber orientation, in shear-out ultimate failure. 

A similar failure behavior is observed near bolt 3, as 
presented in Table 7. At load increment 24, corre- 
sponding to a load of 13,965 lbs, failure occurs in ply 1 
with a 45° fiber orientation, in the net-section failure 
mode, followed by fifteen failures in different plies at 
the same load level until ultimate joint failure. 

Thus, the ultimate joint failure load is computed as 
13,965 lbs near bolt 3. As shown in Table 7, the se- 
quence of ply failure indicates that 90° and ±45° plies 
fail with the net-section failure mode while 0° plies fail 
with the shear-out failure mode. 

Four-Bolt Single-Lap Metal to Composite Joint 

The geometrical parameters for the four-bolt double- 
lap joint shown in Fig. 13 are defined by W= 3.125 in, 
s = 1.25 in, e = 0.9375 in, £ = 2.75 in, and D = 0.3125 
in. 
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Fig. 17 Variation of ply loads through the joint thickness. 
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Table 2 Progressive ply failure in one-bolt single-lap joint. 


Load 

increment 

Applied joint 
load, (lb) 

Ply 

number 

Ply 

orientation, 

(degree) 

Failure mode 

1 

3656 

10 

90 

net-section 

3 

3693 

11 

90 

net-section 

19 

4287 

11 

90 

net-section ultimate 

22 

4373 

10 

90 

net-section ultimate 

34 

4879 

1 

45 

net-section 

37 

4977 

3 

-45 

net-section 

39 

5027 

5 

45 

net-section 

41 

5077 

2 

0 

bearing 

43 

5128 

4 

0 

bearing 

44 

5128 

7 

-45 

net-section 

46 

5179 

6 

0 

bearing 

48 

5231 

8 

0 

bearing 

50 

5283 

9 

0 

bearing 

51 

5283 

14 

-45 

net-section 

52 

5283 

16 

45 

net-section 

53 

5283 

18 

-45 

net-section 

54 

5283 

20 

45 

net-section 

56 

5336 

1 

45 

net-section ultimate 

57 

5336 

3 

-45 

net-section ultimate 

58 

5336 

5 

45 

net-section ultimate 

59 

5336 

7 

-45 

net-section ultimate 

60 

5336 

12 

0 

bearing 

61 

5336 

13 

0 

bearing 

62 

5336 

14 

-45 

net-section ultimate 

63 

5336 

15 

0 

bearing 

64 

5336 

15 

0 

bearing ultimate 

65 

5336 

12 

0 

bearing ultimate 

66 

5336 

13 

0 

bearing ultimate 

67 

5336 

16 

45 

net-section ultimate 


Table 3 Maximum hole enlargement in a three-bolt double-lap joint. 



Bolt 1, (in) 

Bolt 2, (in) 

Bolt 3 (in) 

Aluminum 

plate 

4.068 x 10 -5 

3.053 x 10" 5 

3.174 x 10" 5 

Composite 

plate 

3.933 x 10' 5 

5.723 x 10' 5 

6.090 x 1 0' 5 


Table 4 Spring stiffness values in a three-bolt double-lap joint. 



Bolt 1 
(lb/in) 

Bolt 2 
(lb/in) 

Bolt 3 
(lb/in) 

Aluminum plate 

40,338 

53,518 

57,991 

Plies in composite plate 

45° 

29,071 

22,589 

28,862 

0° 

54,629 

37,740 

39,221 

-45° 

38,564 

24,102 

21,403 

90° 

9,043 

5,832 

6,433 
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Fig. 18 



- 3 - 2-10123 

Ply loads, /*j\ lb) 

Variation of ply loads through the joint thickness near bolt 3 in a three-bolt joint. 



Fig. 19 Variation of bolt/plate displacement through the joint thickness near bolt 3 in a three-bolt joint. 
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Table 5 Progressive ply failure for bolt 1 in the three-bolt double-lap joint. 


Load 

increment 

Bolt load, 
(ib) 

Applied joint 
load, (lb) 

Ply 

number 

Ply 

orientation, 

(degree) 

Failure mode 

1 

3601 

22402 

19 

0 

shear-out 

4 

3673 

22852 

11 

90 

net-section 

5 

3673 

22852 

17 

0 

shear-out 

7 

3710 

23080 

10 

90 

net-section 

8 

3710 

23080 

15 

0 

shear-out 

10 

3747 

23311 

13 

0 

shear-out 

12 

3784 

23544 

12 

0 

shear-out 

14 

3822 

23780 

8 

0 

shear-out 

15 

3822 

23780 

9 

0 

shear-out 

17 

3860 

24017 

2 

0 

shear-out 

18 

3860 

24017 

4 

0 

shear-out 

19 

3860 

24017 

6 

0 

shear-out 

23 

3977 

24745 

10 

90 

net-section ultimate 

24 

3977 

24745 

14 

-45 

net-section 

25 

3977 

24745 

18 

-45 

net-section 

26 

3977 

24745 

18 

-45 

net-section ultimate 

27 

3977 

24745 

3 

-45 

net-section 

28 

3977 

24745 

3 

-45 

net-section ultimate 

29 

3977 

24745 

7 

-45 

net-section 

30 

3977 

24745 

7 

-45 

net-section ultimate 

31 

3977 

24745 

1 

45 

net-section 

32 

3977 

24745 

1 

45 

net-section ultimate 

33 

3977 

24745 

2 

0 

shear-out ultimate 

34 

3977 

24745 

4 

0 

shear-out ultimate 

35 

3977 

24745 

5 

45 

net-section 

36 

3977 

24745 

5 

45 

net-section ultimate 

37 

3977 

24745 

6 

0 

shear-out ultimate 

38 

3977 

24745 

8 

0 

shear-out ultimate 

39 

3977 

24745 

9 

0 

shear-out ultimate 
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Table 6 Progressive ply failure for bolt 2 in the three-bolt double-lap joint. 


Load 

increment 

Bolt load, 
(ib) 

Applied joint 
load, (lb) 

Ply 

number 

Ply 

orientation, 

(degree) 

Failure mode 

1 

2101 

13124 

11 

90 

net-section 

3 

2122 

13255 

10 

90 

net-section 

17 

2415 

15086 

10 

90 

net-section ultimate 

22 

2513 

15698 

11 

90 

net-section ultimate 

26 

2589 

16174 

18 

-45 

net-section 

29 

2641 

16499 

14 

-45 

net-section 

31 

2668 

16664 

20 

45 

net-section 

33 

2694 

16831 

7 

-45 

net-section 

35 

2721 

16999 

3 

-45 

net-section 

36 

2721 

16999 

16 

45 

net-section 

40 

2804 

17514 

1 

45 

net-section 

41 

2804 

17514 

5 

45 

net-section 

46 

2917 

18225 

18 

-45 

net-section ultimate 

47 

2917 

18225 

3 

-45 

net-section ultimate 

48 

2917 

18225 

5 

45 

net-section ultimate 

49 

2917 

18225 

1 

45 

net-section ultimate 

50 

2917 

18225 

7 

-45 

net-section ultimate 

51 

2917 

18225 

14 

-45 

net-section ultimate 

52 

2917 

18225 

16 

45 

net-section ultimate 

53 

2917 

18225 

12 

0 

shear-out 

54 

2917 

18225 

13 

0 

shear-out 

55 

2917 

18225 

15 

0 

shear-out 

56 

2917 

18225 

17 

0 

shear-out 

57 

2917 

18225 

19 

0 

shear-out 

58 

2917 

18225 

20 

45 

net-section ultimate 

59 

2917 

18225 

2 

0 

shear-out 

60 

2917 

18225 

4 

0 

shear-out 

61 

2917 

18225 

6 

0 

shear-out 

62 

2917 

18225 

8 

0 

shear-out 

63 

2917 

18225 

9 

0 

shear-out 

69 

3066 

19155 

2 

0 

shear-out ultimate 
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Table 7 Progressive ply failure for bolt 3 in the three-bolt double-lap joint. 


Load 

increment 

Bolt load, 
(ib) 

Applied joint 
load, (lb) 

Ply 

number 

Ply 

orientation, 

(degree) 

Failure mode 

1 

2302 

12769 

11 

90 

net-section 

3 

2325 

12897 

10 

90 

net-section 

6 

2372 

13156 

19 

0 

shear-out 

8 

2395 

13288 

17 

0 

shear-out 

10 

2419 

13420 

15 

0 

shear-out 

11 

2419 

13420 

20 

45 

net-section 

13 

2443 

13555 

12 

0 

shear-out 

14 

2443 

13555 

13 

0 

shear-out 

16 

2468 

13690 

6 

0 

shear-out 

17 

2468 

13690 

8 

0 

shear-out 

18 

2468 

13690 

9 

0 

shear-out 

19 

2468 

13690 

16 

45 

net-section 

21 

2493 

13827 

2 

0 

shear-out 

22 

2493 

13827 

4 

0 

shear-out 

24 

2517 

13965 

1 

45 

net-section 

25 

2517 

13965 

5 

45 

net-section 

26 

2517 

13965 

5 

45 

net-section ultimate 

27 

2517 

13965 

1 

45 

net-section ultimate 

28 

2517 

13965 

3 

-45 

net-section 

29 

2517 

13965 

7 

-45 

net-section 

30 

2517 

13965 

10 

90 

net-section ultimate 

31 

2517 

13965 

11 

90 

net-section ultimate 

32 

2517 

13965 

14 

-45 

net-section 

33 

2517 

13965 

14 

-45 

net-section ultimate 

34 

2517 

13965 

2 

0 

shear-out ultimate 

35 

2517 

13965 

3 

-45 

net-section ultimate 

36 

2517 

13965 

4 

0 

shear-out ultimate 

37 

2517 

13965 

6 

0 

shear-out ultimate 

38 

2517 

13965 

7 

-45 

net-section ultimate 

39 

2517 

13965 

8 

0 

shear-out ultimate 


An initial joint load of P- 312.5 lbs is applied to 
the composite laminate while the end of the alumi- 
num plate is constrained. Due to the presence of 
symmetry in geometry and loading, only the results 
concerning bolts 1 and 3 are presented. The maxi- 
mum hole enlargement values associated with these 
bolt holes that were computed from the two-dimen- 
sional analysis are presented in Table 8, and the 
spring stiffness values for each ply are in Table 9. 

The through-the-thickness variation of the ply 
loads near bolt number 1 is shown in Fig. 20. The 
corresponding bolt/hole deformations are depicted in 
Fig. 21. As observed in these figures, the most pro- 
nounced deformation occurs in plies located along 
the plate interfaces. 

As presented in Table 10, the initial ply failure 
near bolt 3 occurs at a load level of 14,507 lbs, in ply 
number 1 with a 45° fiber orientation, in the net- 
section failure mode. The failure progresses with the 


±45° and 90° fiber orientations in the net-section 
mode, and further continues with the failure of plies 
with 0° fiber orientation in the bearing mode. Part of 
the laminate near bolt 1 becomes unstable at load 
increment 50, corresponding to a load of 18,790 lbs, 
in ply 7 with a -45° fiber orientation, in the net- 
section failure. At this load level, seven more failures 
occur in the composite laminate before the laminate 
is assumed to ultimately fail at load increment 71, at 
a load level of 21,599 lbs. Near bolt 3, the initial ply 
initial failure occurs at 8,261 lbs, in ply 2 with 0° 
fiber orientation, in the shear-out failure mode as 
presented in Table 11. The failure progresses with 
plies of 0° fiber orientation in shear-out mode. Part of 
the laminate near bolt 3 becomes unstable at load 
increment 29, corresponding to a load of 9,401 lbs, in 
ply 3 with a -45° fiber orientation, in the net-section 
ultimate failure. At this load level, 1 1 more failures 
occur in the composite laminate before the laminate 
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Table 8 Maximum hole enlargement in the four- 
bolt single-lap joint. 



Bolt 1, (in) 

Bolt 3, (in) 

Aluminum 

plate 

0.4144 x IO 4 

0.3169 x 1 O' 4 

Composite 

plate 

0.8532 x IO 4 

1.4378 x 1 O' 4 


Table 9 Spring stiffness values in the four-bolt 
single-lap joint. 



Bolt 1, (lb/in) 

Bolt 3, (lb/in) 

Aluminum plate 

39,527 

63,183 

Plies in 
composite 
plate 

45° 

37,955 

21,018 

0° 

51,811 

38,542 

i 

o 

34,835 

29,017 

o 

O 

OS 

8,747 

6,516 


% 
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.3 
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Aluminum plate 


20 


■: 15 
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30 = 
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20 := 


15 = 
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5 II 


Composite laminate 
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Fig. 20 Variation of ply loads through the joint thickness near bolt 3 in the four-bolt joint. 
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Fig. 21 Variation of bolt/plate displacement through the joint thickness near bolt 3 in the four-bolt joint. 


Table 10 Progressive ply failure for bolt 1 in the four-bolt single-lap joint. 


Load 

increment 

Bolt load, 
(ib) 

Applied joint 
load, (lb) 

Ply 

number 

Ply 

orientation, 

(degree) 

Failure mode 

1 

3296 

14507 

1 

45 

net-section 

6 

3430 

15096 

5 

45 

net-section 

14 

3678 

16185 

16 

45 

net-section 

15 

3678 

16185 

20 

45 

net-section 

20 

3827 

16842 

10 

90 

net-section 

21 

3827 

16842 

11 

90 

net-section 

24 

3904 

17181 

1 

45 

net-section ultimate 

25 

3904 

17181 

5 

45 

net-section ultimate 

26 

3904 

17181 

10 

90 

net-section ultimate 

27 

3904 

17181 

16 

45 

net-section ultimate 

28 

3904 

17181 

11 

90 

net-section ultimate 

29 

3904 

17181 

20 

45 

net-section ultimate 

32 

3984 

17526 

2 

0 

bearing 

34 

4022 

17701 

4 

0 

bearing 

37 

4103 

18057 

6 

0 

bearing 

39 

4144 

18238 

8 

0 

bearing 

40 

4144 

18238 

9 

0 

bearing 

42 

4185 

18420 

3 

-45 

net-section 

43 

4185 

18420 

12 

0 

bearing 

45 

4227 

18604 

13 

0 

bearing 

46 

4227 

18604 

15 

0 

bearing 

47 

4227 

18604 

17 

0 

bearing 

48 

4227 

18604 

19 

0 

bearing 

50 

4270 

18790 

7 

-45 

net-section 

51 

4270 

18790 

7 

-45 

net-section ultimate 

52 

4270 

18790 

3 

-45 

net-section ultimate 

53 

4270 

18790 

14 

-45 

net-section 

54 

4270 

18790 

14 

-45 

net-section ultimate 

55 

4270 

18790 

18 

-45 

net-section 

56 

4270 

18790 

18 

-45 

net-section ultimate 

71 

4908 

21599 

2 

0 

bearing ultimate 
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Table 11 Progressive ply failure for bolt 3 in the four-bolt single-lap joint. 


Load 

increment 

Bolt load, 
(ib) 

Applied joint 
load, (lb) 

Ply 

number 

Ply 

orientation, 

(degree) 

Failure mode 

1 

2260 

8261 

2 

0 

shear-out 

4 

2305 

8427 

4 

0 

shear-out 

6 

2329 

8511 

6 

0 

shear-out 

8 

2352 

8596 

8 

0 

shear-out 

10 

2375 

8682 

9 

0 

shear-out 

12 

2399 

8769 

12 

0 

shear-out 

13 

2399 

8769 

13 

0 

shear-out 

14 

2399 

8769 

15 

0 

shear-out 

16 

2423 

8857 

10 

90 

net-section 

17 

2423 

8857 

17 

0 

shear-out 

18 

2423 

8857 

19 

0 

shear-out 

20 

2447 

8945 

11 

90 

net-section 

21 

2447 

8945 

11 

90 

net-section ultimate 

25 

2521 

9216 

3 

-45 

net-section 

27 

2547 

9308 

7 

-45 

net-section 

29 

2572 

9402 

3 

-45 

net-section ultimate 

30 

2572 

9402 

7 

-45 

net-section ultimate 

31 

2572 

9402 

10 

90 

net-section ultimate 

32 

2572 

9402 

14 

-45 

net-section 

33 

2572 

9402 

14 

-45 

net-section ultimate 

34 

2572 

9402 

1 

45 

net-section 

35 

2572 

9402 

1 

45 

net-section ultimate 

36 

2572 

9402 

2 

0 

shear-out ultimate 

37 

2572 

9402 

4 

0 

shear-out ultimate 

38 

2572 

9402 

5 

45 

net-section 

39 

2572 

9402 

5 

45 

net-section ultimate 

40 

2572 

9402 

6 

0 

shear-out ultimate 

41 

2572 

9402 

8 

0 

shear-out ultimate 


is assumed to ultimately fail at load increment 41. 
Thus, the ultimate joint failure load occurs near bolt 3 
at a load level of 9,402 lbs. 

Conclusions 

In this study, an approach to predict the strength of 
single- and double-lap bolted composites has been 
developed based on the through-the-thickness ply 
loads of the laminate in conjunction with the average 
stress failure criterion. This approach utilizes the 
model of a beam on an elastic foundation to compute 
the corrected ply loads utilizing a two- dimensional 
stress analysis based on the complex potential and 
variational formulation. In the case of a one-bolt 
single-lap aluminum-to-composite joint, the joint 
strength prediction from the present approach is in 
acceptable agreement with the experimental measure- 
ment published previously. This approach proves that 
the ply load distribution in a laminate is significantly 
influenced near the bolt by the bolt bending deforma- 


tions. This distribution is dependent on the plate 
thickness and laminate lay-up, and it is different for 
single- and double-lap bolted joints. 
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Appendix 

In the k th laminate, the strain energy of the / th bolt, which is defined by a uniform cross-section, A (£) , moment 
of inertia, I (0 , and Young’s and shear moduli, E (£) and G (£) , respectively, can be expressed as 


Bf = Y-q< ,0 VV 
T! 2 


(Al) 


where qf represents the vector of nodal deflections and rotations for the i th beam element and K is the number of 
nodes in the bolt discretization. The stiffness matrix for a two-node Timoshenko beam element is given by Ghali and 
Neville 7 as 
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(A2) 


where a ( k lf) = (h% l) ) 2 G {1) A {1) /\2E {l) I {l) c in which c represents the shear correction factor. 

Rearranging the right-hand side of Eq. (Al) such that the matrices are suitable for static condensation of the 
internal nodal rotations leads to 
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where 




= {a::: c 


U) I 


) 

Y k,(K- 1) J 




(A3) 


(A4) 


20 

American Institute of Aeronautics and Astronautics 



b ' . = 


'u ><> I /,( «) I /,<«> I 

_mi_ ___yi 1 

7do "fT?i(j“| Too 1 




A: ,13 


1 # 


A, 33 

oo" 


1 c +c 

| k,22 k,\l 

i 

r b m 

| k, 12 

1 

n 

i 

i 

i 

T 

1 

1 

1 

i 

L 

1 

_L 

i 

j 

1 

J, 

1 

1 

U 

1 J-(( K ~2)£) 

| « ,12 

1 , ((*-2)0 , ((*-1)0 1 
| b k,22 +b k, 11 \ 




b ( = 


U k, 14 


b (U) +b (2£) 

U k , 24 ^ Vt,13 


£ (2 -0 

U k,23 


k,23 


L «*- \)£) 
U k, 12 

~7(~K~\~£~ 

°k,\4 


j foHK-iw 


k, 12 A, 14 , 

T((^-i)“f I 7«F-ijij 

D k,22 I ”*,24 

rr (( ^-i)7rti«^T“[ 

1 _ ''ii 4 ! l | 


& (2/) 

fr,14 


+c 

A, 24 k, 13 


^*,44 ^ U k,33 


b^ 

U k, 34 


*A,44 ~ U k,33 


b^ 

U k, 14 


& (3/) 

U k,34 


' 


U k, 23 


lALK- 2)t) 

U k, 14 

~r([K-2~£~~~f7(K~ 
U k , 24 " l " ‘A, 13 


, ((*-1)*) 
U k, 23 


h«K-\)£) 

U k,34 


i{{K-3)£) i((K- 2)1) 

J k, 44 " + " k,33 

'7({K~2)i) 

*■,34 


l((K-2)£) 


+ b\ 


((*- 1)0 


(A5) 


(A6) 


(A7) 


Since the internal nodes are not subjected to external moments, the first variation of the strain energy with 
respect to the vector q‘^ vanishes, resulting in the moment equilibrium equations as 


b ( ° r A ( / ) +b ,0 q < ' 

k,A $ k k,</></>*-k , 


= 0 


(A8) 


Solving for q^ in the Eq. (A8) and substituting into Eq. (A3), after rearranging the terms, leads to 


-A ( : )T b\ n A\ ( 


where the matrix b[ e) is defined as 
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The strain energy of the translational and rotational spring elements can be expressed as 


s;’=-a: 


(i)T J (?) . (0 


k A, 


(A9) 
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(All) 


in which A‘° is the same as \[ e) with an additional degree of freedom representing the rigid base for the springs, and 
k k ] is defined as 
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where the last column and row appear due to degrees of freedom associated with the rigid base. Since the rigid base 
is fixed, this column and row are eliminated from the matrix in Eq. (A 12). Thus, the total strain energy of the beam 
on an elastic foundation can be written as 


or 


in which k ( k n is defined by 
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(A13a) 

(A13b) 


(A 14) 


In the case of a single-lap bolted joint, the bolt nut is not present in the top plate, k= 1, while the bolt head is not 
present in the bottom plate, k =2; thus, k\ f) n - 0 in the top plate and k ( 2 l) h = 0 in the bottom plate. In the absence of 

both head and nut in the middle plate of a double-lap bolted joint, in addition to the values of k\ l) n - kf ] h - 0 for the 
top and bottom plates, k 2 \ - k { 2 \ - 0 in the middle plate. 
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The total potential energy of the I th bolt can be expressed as 


where 


and 
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(A 15) 


(A 16) 


(A 17) 


for a double-lap joint. The variables K , M , and P define the number of nodes in the beam discretization for each 
plate. 

According to the bilinear behavior of the ply loads shown in Fig. 10, the coefficients of the matrix k[ l) become 
kjfl [ k^) - ak [ l ] , where a - 0.1 ] for initially damaged plies and zero for totally damaged plies. Also, for initially 
damaged plies, the values of (k^j - k^J )^J N are added to the corresponding locations in the vector of externally 
applied loads, F . 

Appling the principle of minimum potential energy and forcing the first variation of the total potential to vanish, 
8n = 0 , leads to the system of equilibrium equations 
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K (/) A (/) = F W (A 19) 

In the finite element formulation of the problem, additional constraint equations are introduced to ensure the 
continuity of the rotational degree of freedom between the adjacent section of the beam 

= &kI for a single lap joint 

(A20) 

i A+1 for a double lap joint 

=</> {e) 

k+M ' K + M + 1 

These constraint conditions are enforced by adding an additional row and column with all but two elements set 
to zero in the matrix K (0 . Two nonzero elements are set to 1 and -1 in the locations corresponding to the nodal 
rotations; a zero value is added in the right-hand side. 
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